Microscopic theory of the magnetoresistance of disordered superconducting films 
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Experiments on disordered superconducting thin films have revealed a magnetoresistance peak of 
several orders of magnitude. Starting from the disordered negative-!/ Hubbard model, we employ an 
ab initio approach that includes thermal fluctuations to calculate the resistance, and fully reproduces 
the experimental phenomenology. Maps of the microscopic current flow and local potential allow 
us to pinpoint the source of the magnetoresistance peak - the conducting weak links change from 
normal on the low-field side of the peak to superconducting on the high-field side. Finally, we 
formulate a simple one-dimensional model to demonstrate how small superconducting regions will 
act as weak links in such a disordered thin film. 
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The interplay of disorder and superconductivity has 
been a subject of great interest and debate for many 
years [TH3]. While the consequences of disorder are rel- 
atively well understood within BCS mean-field theory, 
the way disorder modifies the magnetic field and the 
temperature response of low-dimensional superconduc- 
tors (SCRs), where the loss of superconductivity is due 
to phase fiuctuations [3^ , remains an open question. This 
situation is underscored by recent puzzling experiments 
that demonstrated a huge magnetoresistance (MR) peak 
on the normal side of the SCR-insulator transition in thin 
films [5j , the emergence of a "super-insulating" phase [61 
in the same system, and the smearing of the transition 
at the two-dimensional SC interface formed between two 
insulating oxides [7j . One of the main reasons for our lim- 
ited understanding of these experiments is that there is 
no microscopic theory that incorporates disorder, phase 
fluctuations and magnetic field. In this Letter we utilize a 
recently developed ah initio formalism fS^ to address the 
puzzle of the MR peak. This tool allows us to calculate 
the conductance through a possibly disordered supercon- 
ducting (SC) region, based on a microscopic model, using 
an approximation that takes into account the phase and 
amplitude fluctuations of the SC order parameter. More- 
over, the detailed information provided by this approach 
on the local currents and chemical potentials (for details 
see [S]), in addition to the nature of the current flow 
(electrons or Cooper pairs) illuminates the microscopic 
physics behind the anomalous resistance peak. 

The SC region is described by the negative-?/ Hubbard 
model, a lattice model that includes on-site attraction, 
and may include disorder, orbital and Zeeman magnetic 
fields, and Coulomb repulsive interactions. The Hamil- 
tonian is 

i,cr 

- E U^cl^4iC^lC,t + E Up^'f'^nJ , (1) 

where is the on-site energy, tij the hopping element 
between adjacent sites i and j, and Ui is the onsite two- 



particle attraction. Both \tij \ — t and Ui = U are taken 
to be uniform. An orbital magnetic field can be incor- 
porated into the phases of the hopping elements tij (a 
Zeeman field, manifested as spin-dependent on-site ener- 
gies Eicr, will not be included here). To account for dis- 
order, Ci will be drawn from a Gaussian distribution of 
width W. The local density is fii = '^icr'^icj ^-i^d the 
long-range repulsion is given by the screened Coulomb 
interaction, J7j^°"' — J/c exp[— rij/AJa/ri^ , with Tij the 
distance between sites i and j, and A = 2a the screening 
length, both measured in units of the lattice constant 
a (for computational efficiency we cut off this interac- 
tion after four sites). Unlike, for example, the disordered 
XY model, the negative-!/ Hubbard model can lead to a 
BCS transition, a BKT transition [4 , or to a percolation 
transition and thus this choice will not limit a priori 
the underlying physical processes. In order to compare 
the results of the model to experiments, one needs to 
compare the relevant length scales in particular the SC 
coherence length and the mean-free path (or localization 
length) . In the calculations presented below we keep the 
former (which is governed practically by U/t) fixed, and 
tune the latter by changing the disorder W/t. 

In Ref.[8] we have developed an exact formula for the 
current through such a SC region, sandwiched between 
metallic leads. Furthermore, the resulting expression for 
the current can be evaluated using a Monte Carlo ap- 
proach that treats the fermionic degrees of freedom ex- 
actly, takes into account the thermal fluctuations of the 
amplitude and phase of the SC order parameter, but 
neglects its quantum fluctuations. It has been demon- 
strated [8] that such an approach describes quantitatively 
various transport processes in the vicinity of the BKT 
transition, which are crucial to describe the physics of 
low-dimensional superconductors . 

In the following we first recover the main experimen- 
tal results and study how the microscopic variables af- 
fect the resistance. We then construct a detailed map of 
the current flow through the system to uncover the mi- 
croscopic mechanism that drives the MR peak. First, in 
Fig. [IJa) we depict the resistance at several temperatures 
for a sample of size Lx = 21a {x is the direction of the 
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FIG. 1: (Color online) (a,b) The MR with an applied mag- 
netic field (a) perpendicular to the surface and (b) parallel to 
the surface. The blue curve (highest at large fields) is at low 
temperature, yellow intermediate and red high temperature. 
The arrows highlight the critical field Be and the maximum 
resistance field Bm- The points (i - vi) denote the magnetic 
fields where the current maps of Fig. [S] were evaluated. The 
perpendicular MR for four different amplitudes of disorder is 
shown in (c), and with Coulomb repulsion in (d). 



applied potential), Ly = 11a, and = 2a. The interac- 
tion strength is U = 2t, the onsite disorder has width 
W — 0.5t, the average density is at 41% filling, and 
initially there are no long-range repulsive interactions, 
Uc = 0- Throughout we set fee = 1- At zero magnetic 
field, the resistance of the SC state is solely the lead-SCR 
contact resistance. Increasing the magnetic field perpen- 
dicular to the surface drives a SCR-insulator transition, 
around Be — 0.04, in units of (electron) flux quantum 
per square [hc/ea^], manifested in the reversal of the tem- 
perature dependence of the resistance. Above the transi- 
tion the system displays, in agreement with experiment 
[5], a peak resistance (at the field Bm — 0.07) that is 
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FIG. 2: (Color online) (a) The temperature dependence of 
the resistance at three values of the magnetic field, (b) The 
activation energy Ta with magnetic field, (c) The activation 
energy Ta (pluses, red), superconducting transition temper- 
ature Tc (crosses, green), and ratio of mean free path £ to 
coherence length ^ (squares, blue) with disorder, (d) Ta vs 
Tc and screened Coulomb interaction Uc for different disorder 
levels. 



over a hundred times the normal state resistance. This 
feature consistently emerges for several different statis- 
tical configurations of disorder, with the peak resistance 
varying by a factor of approximately ten. On changing 
the direction of the (orbital) magnetic field from normal 
to lying in the plane of the sample (Fig. [ijb)), both Be 
and Bm shift to higher fields, with a slight decrease in 
the magnitude of the MR peak, in complete agreement 
with recent experiments ^lOJ. 

We now explore how the MR varies with the parame- 
ters of the system. Firstly, in Fig. [ijc) we focus on one 
specific realization of disorder but increase its amplitude. 
At low disorder, W — O.lt, the magnetic field suppresses 
the SC state, but there is no MR peak. On increasing 
the magnitude of the disorder to W = 0.2t a peak in the 
resistive curve emerges at B ^ 0.09. (Interestingly, the 
peak emerges approximately when the mean free path £ 
becomes smaller than the coherence length ^, Fig. [2jc).) 
As disorder is further increased both the SC transition 
and MR peak shift to lower magnetic fields [TT] . In agree- 
ment with experiment, the MR peak persists even when 
the SC phase is suppressed at zero field {W — 0.5t and 
W = t), until the peak is extinguished at large enough 
disorder (W = 2t). Fig. [l|^d) demonstrates that the MR 
peak becomes more dramatic with increasing Coulomb 
repulsion, as the MR peak increases by two orders of mag- 
nitude. The enhancement of the MR peak by both dis- 
order and Coulomb repulsion points towards a Coulomb 
blockade mechanism. 

Fig. depicts the temperature dependence of the 
resistance for several values of magnetic field, which 
demonstrate an activated behavior, R = i?o exp(rA/T). 



3 




(iv) B = 0.09fec/ea^ T = 0.04t (v) B = 0.08/tc/ea^ T = 0.02t (vi) B = 0.08fee/ea^ T = 0.08t 




FIG. 3: (Color online) Current maps at six values of the normal magnetic field shown in Fig.jTJa). The net current flow is from 
left to right. The cyan darts show supercurrent and magenta pointers normal current direction and magnitude; their magnitude 
is renormalized in each map according to the maximum local current flow. The white lines show equi-chemical potential lines 
across the sample (ten lines correspond to the full potential drop between the two leads), and the background shading depicts 
the magnitude of the local superconducting order parameter ranging from strong (orange) to zero (black). The labels (1) to 
(3) highlight three points of special interest. 



The activated behavior, and the peak of the activation 
temperature Ta at the field Bm (Fig. [2jb)) agrees with 
experimental observations [S]. Ta increases with disor- 
der (Fig.|2jc)) and with repulsive interactions (Fig.[2]^d)), 
consistent with the enhancement of the MR peak by these 
parameters (Fig.[l|c) and (d)). The critical temperature 
Tc at which SC correlations are suppressed at zero-field 
on the other hand falls with disorder. In Fig. [2jd) this 
gives us an inverse linear dependence of Ta on Tc, similar 
to that seen in experiment [S]. 

So far we have demonstrated that the calculation re- 
produces the phenomenology of the experimental obser- 
vations and hinted that activated transport through a 
blockade region drives the rise of the resistance. We now 
take advantage of our ability to probe the evolution of 
the local currents and chemical potentials with applied 
magnetic field (Fig. [S] and comprehensively in Fig. S2) 
to pinpoint the physical processes behind this behavior. 
Below the critical field. Be — 0.04, Fig. S2 shows that 
no voltage drops across the sample, though the current 
flow, consisting of Cooper pairs, becomes nonuniform. 
As the field increases, {B = 0.06, Fig. |3ji)), there is no 
longer long-range SC coherence across the system, and 
current crossing the system changes its nature between 
SC and normal. Above this field the current has to tra- 
verse normal areas of the sample, and as their density 
increases with magnetic field, the resistance rises. As the 
field increases further {B = 0.07, Panel (ii)), the disor- 
der induces specific channels of transport - the current 
starts as SC, changes into normal current around point 
(1) and then splits towards points (2) and (3) where it 
reverts to Cooper pairs, only to change back into normal 
current as it enters the right-hand lead. At this magnetic 
field, near the maximal value of the resistance, there are 



approximately equal contributions to the current from 
electrons and from Cooper pairs. With increasing field 
{B = 0.08, Panel (in)) the SC areas (2) and (3) shrink 
and become the main source of voltage drop and resis- 
tance in the sample. A larger field {B = 0.09, Panel (iv)) 
suppresses the SC correlations, lowering the resistance of 
the weak links and the overall the resistance. This de- 
crease is, in fact, an interplay of two phenomena ~ larger 
SC areas do not serve as weak links (see point (1) in 
Fig. |3jiii)), but as they shrink with increasing magnetic 
field the resistance associated with them increases sig- 
nificantly (see point (1) in Fig.[3];iv)). At the same time 
small SC areas, that gave a significant contribution to the 
resistance (point (3) in these two panels) become normal 
and the overall resistance decreases. A much larger mag- 
netic field (Fig. S2) suppresses superconductivity almost 
completely, giving rise to a uniform drop of voltage across 
the sample and a significantly lower resistance. Similar 
behavior has been obtained in the presence of repulsive 
Coulomb interactions (not shown). 

To explore the origin of the Arrhenius behavior pan- 
els (v) and (vi) depict the current flow a,t B = 0.08 at 
lower and higher temperatures than Panel (iii) . At point 
(1) we see direct evidence of the activated transport: on 
lowering the temperature (Panel (v)) more potential is 
dropped across the weak SC link, whereas with increas- 
ing temperature (Panel (vi)) the potential dropped across 
this weak link falls. The effect of the weak link is so 
profound that in the hotter panel (vi) despite the reduc- 
ing SC order more supercurrent flows as the electrons no 
longer skirt around the SC island as they did in Panel 
(v). Interestingly at other places in the sample we see 
the more conventional effect of temperature suppressing 
superconductivity, for example at point (3) the increas- 
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FIG. 4: (Color online) The inset of (a) shows the ID model of 
a normal state (blue spheres) conducting channel with weak 
disorder containing a superconducting grain (red prisms) of 
varying length Lsc- In (a) to (c) the red dotted curve is for the 
same set of standard parameters, and the green solid curves 
correspond to (a) higher dot order parameter A = 0.9t, (b) 
higher temperature T — O.It, and (c) with Coulomb repulsion 
Uc = 0.02. 

ing temperature from Panel (v) to Panel (ill) reduces the 
SC order parameter and the supercurrent. 

Why would small SC areas serve as weak links and con- 
tribute to the increasing resistance? To shed light on this 
question we studied a toy-model, consisting of a SC is- 
land with finite order parameter A = 0.3t of length Lsc, 
embedded in disordered normal chain with A = 0, equal 
nearest neighbor tunneling, and a temperature T = 0.05i, 
shown in the inset of Fig. [4][a) . The main figure depicts 
the resistance of the chain as a function of Lsc, with 
the total chain length fixed. As the SC region becomes 
shorter, there is a mild increase in the resistance due the 
increasing length of the resistive normal regions. How- 
ever, when the SC segment is shorter than 4a (similar 
in size to the SC islands around the resistance peak in 
Fig. [sj^iii)), there is a sudden increase in the resistance 
of the sample. The resistance peak is enhanced by rais- 
ing the SC order to A = 0.9t, (Fig. l4]^a)), or increas- 
ing temperature to T = O.lt (Fig. |4]jb)), demonstrating 
an activated transport process. This is enhanced on ap- 
proaching the Anderson limit, where the single particle 



level-spacing in the SC segment becomes larger than A 
[121 . Further corroboration of this picture comes from 
Fig. ^a) where the resistance peak emerges at shorter 
Lsc for higher A. Introducing repulsive interactions 
Uc = O.Obt (Fig.|4];c)) has two major effects - there is a 
larger increase in the resistance and the effect starts at 
much larger SC segment size, because now the SC gap has 
to be compared to the Coulomb blockade energy rather 
than the single-particle spacing. Note that this happens 
even though this segment is not weakly coupled to the 
rest of the chain. This is consistent with the description 
of the MR peak in the two-dimensional systems, where 
as the SC islands shrink the level spacing, or, more re- 
alistically the Coulomb blockade energy becomes of the 
order of the SC gap, the islands become the weak links 
and drive up the resistance of the sample. 

To conclude, we have demonstrated the emergence of 
a magnetoresistance peak, starting from a microscopic 
negative-?/ Hubbard model. The resistance peak is 
driven by the competition between Cooper pair and elec- 
tron transport, where on the weak field side of the peak 
the normal regions serve as the weak links across the sam- 
ple, while this role is played by the SC islands on the high 
field side of the peak, where Coulomb repulsion proba- 
bly plays a dominant role. These ideas can be further 
tested by checking the effect of screening by a parallel 
metallic gate on the peak characteristics. Our picture 
has some similarities to the heuristic model presented in 
Ref. [13] , though there the resistance on the high- field side 
of the peak was completely associated with the normal 
electrons, while it is different than the phenomenologi- 
cal approach of Refs. [THfT5] that discussed the interplay 
between quasi-particle and Cooper pair transport. 
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FIG. 5: The variation of resistance with an applied magnetic 
field normal to the surface at T = 0.04t. Ro — h/e^ is the 
resistance quantum. The points (i - xi) denote the magnetic 
fields where the current maps of Fig. |6] were evaluated, Be 
represents the magnetic field where the temperature depen- 
dence of the resistance changes sign and Bm the magnetic 
field at the resistance maximum. 



I. SUPPLEMENTARY MATERIAL FOR 
"MICROSCOPIC THEORY OF THE 
MAGNETORESISTANCE OF DISORDERED 
SUPERCONDUCTING FILMS" 

In this supplement we present a comprehensive series of 
current maps to fully expose the microscopic mechanism 
driving the breakdown of the superconductor (SC) with 
applied magnetic field shown in Fig. [5] The maps shown 
in Fig. [6] cover the entire scope of the SC-insulator tran- 
sition: maps (i - ii) capture the slight rise of resistance 
as the superconducting order parameter is reduced but 
the superconductor remains phase coherent. Maps (iii - 
vii) outline the loss of phase coherence and formation of 
weak links through the sample, which drives a dramatic 
increase in the sample resistance. Finally, in maps (viii - 
xi) we witness the last remnants of superconductivity be- 
ing destroyed and the system entering the normal phase 
with the resistance plateauing out. 

We first study the maps (i - ii) that show the current 
flow when the superconductor is still phase coherent and 
the resistance is low. The first map (i), corresponds to 
B = 0.01, below the critical field Be (the unit of the 
magnetic field is the electron quantum flux, hc/e, per 
lattice unit area a^). The supercurrent flow is uniform 
and no voltage is dropped across the sample, but instead 
the potential is dropped across the lead-superconductor 
junction. As the field increases to B = 0.04, (map (n)). 



vortices penetrate the sample at points of strong disorder. 
The supercurrent is no longer uniform as it encircles the 
vortices, but the sample remains phase coherent and still 
no potential is dropped across the sample. 

Above the critical field Bq, the system loses long-range 
SC coherence and the current changes its nature between 
SC and normal as it crosses the system. As the field in- 
creases to B = 0.05 (map (iii)) we see the first signs of 
decay of the order parameter and the current flow be- 
comes normal. With the loss of phase coherence across 
the system, for the first time potential is dropped across 
the sample, though around half of the total potential dif- 
ference is still dropped across the lead-superconductor 
barriers. Around this field the temperature dependence 
of the resistance changes sign from superconducting to 
insulating. As the field increases to B = 0.06, (map 

(iv) ), normal regions start to form in the sample, and 
the potential is now dropped mainly across these weak 
links. As the field increases further to B — 0.07, (map 

(v) ), the disorder induces specific channels of transport - 
the current map in Fig. |6] reveals that the current starts 
as SC, changes into normal current around point (1) and 
then splits towards points (2) and (3) where it reverts 
to Cooper pairs, only to change back into normal cur- 
rent as it connects to the right lead. At this magnetic 
field Bm, the maximum in the resistance, there are ap- 
proximately equal contributions to the current from elec- 
trons and from Cooper pairs. On increasing the field to 
B = 0.08, (map (iv)), the SC areas (2) and (3) shrink and 
become the main source of voltage drop and resistance 
in the sample. A larger field B = 0.09, (map (v)), sup- 
presses the SC correlations, making area (3) normal, low- 
ering the resistance of that weak link, while at the same 
time, the SC area (1) shrinks, giving rise to a higher local 
resistance. Overall the resistance of the sample decreases 
due to the disappearance of some of the SC island weak 
links. 

Finally, above B — 0.10 (maps (viii) and (ix)) the 
SC islands shrink to such extent that the current fiows 
predominantly around them, through the normal areas. 
With a further increase in the magnetic field to i? = 0.12 
(map (x)) only normal current flows through the channels 
that avoid the residual weak pairing. However, a further 
increase in magnetic fleld to _B = 0.15, (map (vi)) sup- 
presses superconductivity almost completely, giving rise 
to a uniform drop of potential across the sample. 
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(iv) B = 0.06fec/eo^ 
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FIG. 6: Current maps at eleven values of the magnetic field shown in Fig. [5] The net flow of current is from left to right. 
As illustrated in the key, the cyan darts show supercurrent and magenta pointers normal current direction and magnitude; 
their magnitude is renormalized in each map according to the maximum local current flow. The white lines show equi-chemical 
potential lines across the sample (ten lines correspond to the full potential between the two leads), and background shading 
depicts the magnitude of the local superconducting order parameter ranging from strong (orange) to zero (black). 



